Nonequilibrium Cotunneling through a Three-Level Quantum Dot 
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We calculate the nonlinear cotunneling conductance through a quantum dot with 3 electrons 
occupying the three highest lying energy levels. Starting from a 3-orbital Anderson model, we apply 
a generalized Schrieffer- Wolff transformation to derive an effective Kondo model for the system. 
Within this model we calculate the nonequilibrium occupation numbers and the corresponding 
cotunneling current to leading order in the exchange couplings. We identify the inelastic cotunneling 
thresholds and their splittings with applied magnetic field, and make a qualitative comparison to 
recent experimental data on carbon nanotube and InAs quantum-wire quantum dots. Further 
predictions of the model such as cascade resonances and a magnetic-field dependence of the orbital 
level splitting are not yet observed but within reach of recent experimental work on carbon nanotube 
and InAs nanowire quantum dots. 
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The Kondo effect has been observed in a number of 
different quantum dot (QD), and single molecule de- 
vicesii^i^i^. The effect is manifested as a sharp con- 
ductance peak at zero bias-voltage, developing when the 
temperature is lowered beyond the characteristic Kondo 
temperature. It relies on a spin-degenerate ground-state 
on the quantum dot, which gives rise to logarithmically 
singular spin-flip scattering of the conduction electrons 
traversing the dot. Meanwhile, for quantum dots with 
sufficiently small level spacings or slightly broken de- 
generacies, the Kondo-peak at zero bias voltage can be 
flanked by two or more satellite steps or even peaks 
in the nonlinear conductance. Such inelastic cotunnel- 
ing features are often seen in both GaAs^, Carbon nan- 
otube'''^'^ (CNT) and InAs-wire'^'^" quantum dots, but in 
most cases they are masked by charge-excitations which 
can be nearby in energy and for this reason they have 
not received much attention. 

Single-molecule transistors, on the other hand, exhibit 
a much larger charging energy, Ec 100 meV instead 
of 5 meV, say, for a typical quantum dot). At the same 
time, these molecular systems often display a number of 
degeneracies which are weakly broken once the molecule 
is contacted by source, and drain electrodes, thus induc- 
ing splittings of the order of a few meV. Most recently, 
this was seen in Refs. [TTHl2l . where junctions holding an 
0PV5, or a Ceo molecule, respectively, showed a very 
clear singlet-triplet splitting on the scale of 1 meV to- 
gether with a charging energy of the order of 100 meV. 
This is a very convenient separation of energy scales 
which raises the experimental resolution of inelastic co- 
tunneling phenomena to new standards. Nevertheless, 
given the simpler level-structure of most conventional 
QD-devices, it is desirable to revisit and understand the 



details of interorbital transitions better in these systems. 
This is what we set out to do in the present paper. 

In the following, we study the case of a quantum dot 
occupied by an odd number of electrons, featuring a spin- 
doublet groundstate, giving rise to a zero-bias Kondo 
peak, but with additional orbitals/levels leading to flank- 
ing inelastic cotunneling steps or peaks. The basic 3- 
orbital Anderson model is illustrated in Fig. [1] and will 
be shown to host a variety of different I-V-characteristics 
depending on the relative magnitudes of the 6 different 
tunneling- amplitudes. This 3-orbital model was brieffy 
discussed by some of the present authors in Ref . (cf. in- 
set in Fig. 4b), where it was invoked to explain the rel- 
atively sharp peaks at finite bias flanking a zero-bias 
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FIG. 1: Illustration of the setup. We consider the three 
highest-lying levels of a quantum dot, separated by energies 5, 
and 5' . Both orbitals are connected to source (left lead) and 
drain (right lead) via 6 different tunnel couplings tcm- The 
dot is in the Coulomb blockade regime with the three levels 
adjusted by a gate to accommodate exactly 3 electrons. 
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FIG. 2: (color online) Differential conductance, G, as a func- 
tion of bias voltage, Vsd, for an InAs-wire based quantum dot 
at T = 0.3 K. The data were taken at magnetic fields (per- 
pendicular to the wire) 5 = (thick), 0.1 (dotted), 0.9 
T (red) and the curves were offset by 0.008 /h for clarity. 
The data were taken for an odd occupied Coulomb diamond 
at gate voltage Vg = —2.35 V"". 



Kondo-effect observed in an InAs quantum-wire dot^. 
The measured nonlinear conductance curves for varying 
applied magnetic fields are shown in Fig. [2] This exper- 
iment constitutes one of the rare cases where such side- 
peaks could actually be resolved. As one further exam- 
ple, Fig. [3] shows similar data recorded on a single- walled 
carbon nanotube (CNT) quantum dot (QD). Both sets 
of measurements show well-defined peaks which split into 
weak thresholds on applying a magnetic field. Notice the 
different field-strengths needed in the two experiments, 
reflecting the roughly 4 times larger spin-orbit enhanced 
g-factor in InAs as compared to the CNT. 

The bulk of this paper deals with leading order 
nonequilibrium cotunncling for the 3-orbital Anderson 
model in the Kondo-regime. First we derive an effective 
cotunneling, or Kondo model for a system with three 
electrons distributed on the three orbitals. From this ef- 
fective low-energy model we then proceed to calculate 
the I-V characteristics to leading (second) order in the 
cotunneling amplitude. In section III we discuss some of 
the salient transport features of this system and in sec- 
tion IV we discuss the CNT and the InAs data shown in 
Figs. [2] and [3l In this context we point out additional 
properties of such samples which can be read off the dl- 
dV characteristic. 



I. EFFECTIVE LOW-ENERGY KONDO MODEL 

The 3-orbital Anderson model corresponding to the 
setup in Fig. [T] we write as 

H — Hicad + Hdot + Htun ■ (1) 
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FIG. 3: (color online) Differential conductance, dl/dV, as a 
function of bias voltage, V , for a carbon nanotube quantum 
dot at T = 0.08 K. The data were taken at magnetic fields 
(perpendicular to the tube) B = Q (thick), 0.1 (dotted), 1 
(thin), 2, 3, 9, 10 T (red) and the curves were offset by 
0.008 e^/h for clarity. The data were taken for an odd oc- 
cupied Coulomb diamond at gate voltage Vg = —4.96 V— . 
(Note that at finite magnetic fields features are broadened 
due to noise induced by the magnet power supply). 

The leads are described by the non-interacting Hamilto- 
nian 

i?lcad= X! ('^k - Aia)4kaCak<T (2) 
a=L.R 
k(T 

where a — L, R labels the leads, assumed to be in equi- 
librium at chemical potentials Ha- The operator cj^j.^ 
creates an electron in lead a of momentum k and spin a. 
A simple constant interaction model is used to describe 
the quantum dot itself in terms of three non-degenerate 
orbitals with one common Coulomb repulsion U{'^ Ec)' 

-ff dot = ^ ^ £na- fna f na- 
na 

+ luJ2Y.fnafna{fL'f,na' ~ l) . (3) 

na ma' 

where n,m — 1,2,3 label the orbitals on the dot and 
creates an electron in orbital n with spin a and with 
energy Cna = f-n - ^(^gfJ-sB where €2 = ei + S, and £3 = 
€2 + 6', with level-splittings denoted by 5 and S' . Finally, 
the tunneling Hamiltonian, 

Htun — ^ ^ {^anC^aaftia ^" h.C.^ , 

a=L.R 
na 

describes the coupling of the leads to the dot via six in- 
dependent tunneling amplitudes tam see illustration in 
Fig. [TJ For later convenience, we have introduced the 
local conduction electron operators cj,^ = Sk'-ako-- 

We now restrict our attention to the Kondo-regime in 
which all three orbitals are sufficiently narrow compared 
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FIG. 4: Illustration of the three lowest lying energy states 
that are taken into account in the effective cotunneling model 
of the 3-orbital Anderson model. From left to right; spin, 
hole and particle states. 



to the charging energy to effectively suppress all charge- 
fluctuations, leaving the dot with a well defined occupa- 
tion of 3 electrons in three orbitals. Since the experimen- 
tal data in Figs. [2] and [3] do not resolve higher lying peaks 
corresponding to excited states with two electrons in or- 
bital 3, we shall simplify our model further by also omit- 
ting these higher lying 3-particle states. In the experi- 
ments, the cotunneling thresholds for these higher lying 
states must be comparable to the charging-energy and 
are therefore masked by charge-fluctuations. As men- 
tioned earlier, this is a typical problem of the relatively 
large quantum dots (compared to single-molecule junc- 
tions) which provide a rather poor separation of energy- 
scales. Altogether, we are now left with an effective 
low-energy Hilbert-space spanned by the six lowest ly- 
ing three-electron states: 



(4) 



and the 4-particle states, that is: 

E^e- =2ei + U> 
E^,- = 2ei + 2£2 + 6U> E^, 



(6) 



In the following, we shall choose £2 — —5/2 U, which 
places the system at the particle-hole symmetric point 
where E^^- - E^^- = E2e- ~ E^^- = U/2, with E^^- = 
Eg^ as the ground state. 

In order to eliminate charge-fluctuations from the 3- 
orbital Anderson model ([1]), we employ a generalized 
Schrieffer- Wolff transformation^^ which serves to elimi- 
nate the tunncling-term, -fftun, and retains only terms of 
second order in the tunneling-amplitudes in the form of 
an effective cotunneling, or Kondo Hamiltonian: 



HcS = i?lcad + ^ Eaa\a<j){aa\ + Hint- 



(7) 



The effective interaction now takes the form of a 
spin/orbital exchange-term: 



Ejab 

af3.ab 
a' a 
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cL'C;3., (8) 



in terms of the vector of Pauli matrices, x, the spin- 
operator for the quantum-dot 



(9) 



and the potential scattering term Pah = P^'^^ + Pal which 
consists of an inelastic scattering which involves a change 
in the orbital state 



^ab'^ —-^{^a.s^b,p + Sa.p^b,s — ^a.shji — Sa,hSb,s) 



J2\aa'){ba'\ 



(10) 



We label these three-body states by indices a,b = 
{s,h,p} (i.e. {spin, hole, particle}), together with the 
spin-index a ~ {Ti J.}- The notation is obvious from the 
illustration of the states in Fig. [H The energies of these 
states are 



Esa - 2ei + £2 + 3C/ - -ag^iBB, 

Eha = El + 2£2 +3U - ^crgfigB, 
1 



(5) 



E, 



pa 



2ei 



3U - -cfq^ibB. 



In the case of equidistant energy levels Eha = Ep^ and 
the two excited states, corresponding to respectively an 
electron moved to orbital 3 or a hole moved to orbital 
1, are degenerate. The energy £2 is set by tuning the 
voltage on a gate and here we choose it such that the 
3-particle state will be lower in energy than both the 2, 



and the elastic scattering which occurs via empty levels, 
i.e. level 3 for \s) and \h) and level 2 for \p), 



Pal —-T^a.b ^ {Sd,sSa,p + Sa,hSa,p + Sa,pSa,s) 



X \aa'){aa'\ 



(11) 



To lowest order a transition between the excited states 
\p) and \h) is not possible. 

A constant energy offset arising in the Schrieffer- Wolff 
transformation is neglected and so are further potential 
scattering terms. These lead only to a constant offset in 
the differential conductance of the order oi S/U"^ , which 
can be neglected in our qualitative study. Details of 
the calculation, including the expression for the coupling 
functions are given in appendix El 

Note that when calculating the effective cotunneling 
amplitudes in the appendix, we retain the differences 
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in energy-denominators of the different amplitudes, i.e. 
at this stage we do not make the approximation that 
U ^ (5'^'-'. The according differences in amphtudes ex- 
presses the broken orbital symmetry of the model, even 
with all six tunneling-amplitudes being equal. Neverthe- 
less, in order for the Schrieffer- Wolff transformation to 
be meaningful, we must demand these differences to be 
small, and in all our numerical calculations we therefore 
choose a large charging energy {U ^ 1005) which effec- 
tively makes all energy denominators equal. In other 
words, retaining S^'"^ in the denominators of the cotunnel- 
ing amplitudes would require a more careful treatment of 
charge-fluctuations . 



II. NONEQUILIBRIUM PERTURBATION 
THEORY 



When the bias voltage is large enough to populate the 
excited states on the dot, these will no longer be ther- 
mally occupiedi^ii^. This effect is incorporated in the 
Keldysh component Dyson equation expressed in terms 
of nonequilibrium Green functionsi^iil with self-energies 
calculated to leading (second) order in the effective co- 
tunneling amplitudes or exchange couplings vpJ^p ■€i 1. 
Following the approach taken in Ref. [l^, we employ a 
pseudo-fermion representation for the dot-states, with 
operators defined by da^|0) = \aa) and (0|d^^ = {acr], 
and subject to the constraint Q = J2a cr '^aa^aa = 1- 
The constraint is enforced with the aid of a Lagrange 
multiplier A, included as an additional term, Hx = 
^J^aa'^aa-'^aaJ ^hc Hamiltonian^* . The exact projec- 
tion to the physical Hilbert space is effected by taking the 
limit X oo. We shall need the non-equilibrium Green's 
bmctions: 

Sab.aa'{T,T') - -i (Tc, (d,, (r)^^, (t')) ) (12) 

Gc;3,..'(r,T') = -i(TcK(c„,(r)ct^,(r'))) (13) 

with Tcf; being the time ordering operator along the 
Keldysh contour. Calligraphic letters denote pseudo- 
fermion, italic letters conduction electron Green's func- 
tions. 

The retarded and advanced Green's functions can be 
calculated like in the equilibrium case directly from the 
Dyson equation. The pseudo-fermion spectral function is 
obtained from the imaginary part of the retarded Green's 
function, 



-2 Im 



1 



aa,a 



(14) 



and takes the approximate form of a Lorentzian at the 
resonance frequency coaa- — £aa + A and of width Taa.a — 
—2 Im[S^jj g.]. Since the spectral function appears in later 
evaluations only in convolution with functions that vary 
on a larger energy scale than Taa.a , we approximate it by 



a simple delta-function: 



(15) 



The lesser function is found from the quantum Boltz- 
mann equation (QBE), or generalized Kadanoff-Baym 
equation: 



(16) 



Within the delta-function (quasiparticle) approximation 
for the spectral function, this equation can be solved us- 
ing the following ansatz: 



SaaA^) = iriaaAaai^^), 



(17) 



through which the QBE takes the form of a simple rate- 
equation. 

The conduction electrons are assumed to remain in 
thermal equilibrium and are therefore characterized sim- 
ply by their respective chemical potentials together with 
a simple flat-band approximation for the momentum- 
summed spectral functions, or local conduction electron 
density of states (DOS) at the contact. 



A{uj) = 2iivfO{D - \lo\), 



(18) 



in terms of the DOS at the Fermi surface, vp, and half 
bandwidth D. 



A. Nonequilibrium occupation numbers 



In writing the QBE in Eq. p6|) . we have tacitly as- 
sumed the pseudo-fermion self-energies to be diagonal in 
both spin, and orbital indices. Neglecting spin-orbit in- 
teractions spin is a conserved quantum number and the 
self-energy will therefore be diagonal in this index. This 
does not hold for the orbital quantum number a and off- 
diagonal terms, Safc.ir: can arise. Nevertheless, in the 
regime studied here for which the level splittings (5, 5' are 
assumed to be much larger than level broadenings. Fa, 
energy conservation suppresses such off-diagonal terms 
and the self energy can safely be assumed to be diagonal. 
For an example in which off-diagonal contributions re- 
main important we refer the interested reader to Ref. [l9l . 

Thus neglecting off-diagonal self-energies and using 
Eq. ini), the QBE in Eq. takes the following form: 



(Watr) 



Taai^-^aa) 



(19) 



The self energy S^^. itself depends on the occupation 
numbers and this therefore constitutes a set of six coupled 
equations for six unknown occupation numbers. These 
equations are underdetermined and should therefore be 
solved together with the constraint Q = 'Yliaa"'^"-'^ ~ ^■ 
More details of the actual calculation of pseudo-fermion 
self-energies are given in appendix |B] 
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B. Cotunneling current 



We obtain the current operator directly from the time- 
derivative of the density operator at the contact in the 
left lead, say (see Ref. [la and references therein) . With 
nL = E^^ctClct' one finds 



duL ie 



J2 l^^"^^"^' 

a,ab,(j'a 



lab t „ jab t „ 



PabSa'.a (^JLR^^^La^Ra- ~ Jrlc\ 



(20) 
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In terms of the correlation functions, 



Dlf{r,r') = -i{Tc^ Ja1S,,(r)ix.vcL(r')c,,(r')) 

ab.a' a 
ab 

Dfkir, r') = Df^\T, r') + Z?r^(r, r'), (21) 

the expectation value of the current is given in lowest 
order in the couplings J£^ by 



47re 



ah,aa' 
1 a' a"^ (j(j' 



(^ 



c={p,/i} 

X (^naCTi^(£h<T' - Eqct - eV) - nba'Y{eaa- - £bcr' + cV) 
d 

(22) 



where Y{x) = x nsix) and nsix) = l/(exp[a::/2r] — 1) 
denotes the Bose distribution. The sum goes over all pos- 
sible states a,b — {s, h,p}, which are weighted according 
to the nonequilibrium occupation numbers riba'- 



III. DISCUSSION OF RESULTS 

Tuning the six different tunneling amplitudes, the two 
different level-spacings and external magnetic field al- 
lows for a large variety of cotunneling I-V characteristics. 
With an eye towards the experiments mentioned in the 
introduction we shall discuss a few of the salient features. 
Since we restrict our calculations to leading order pertur- 
bation theory, we do not address the interesting question 
of Kondo-correlations in neither the zero-bias peak, nor 
any of the finite-bias conductance peaks. 
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FIG. 5: (color online) Occupation numbers of the spin, 
'hole and particle states in the 3-level Anderso n mod el for 
symmetric coupling to the leads, tan = Q.l\/U /uf ■ For 
a clear separation of energy scales the temperature is cho- 
sen small T/5 = 0.001 and the levels are clearly separated 
5' /5 — 2. The dashed lines indicate the transition energies 
{±gfiBB, 5, 5 ignsB, 5' ,5' liz gfiBB} in {blue, black, red} for 
a finite magnetic field of B/S = 0.25. 



A. Occupation numbers 

Results for the occupation numbers are shown in Fig.[5l 
For convenience we perform the calculations with a small 
thermal smearing. Nevertheless, the temperature is cho- 
sen far smaller than any other energy scale in the problem 
(T = 0.001(5 throughout) and all plots can be thought of 
as corresponding to T = 0. For finite magnetic field and 
zero bias- voltage the ground state, \s t), is occupied with 
probability one. As illustrated in Fig. O the other states 
become populated for larger voltages and in the limit 
1/ — > oo all states will be equally occupied. Inter-orbital 
transitions play a role as soon as the voltage becomes 
larger than S or S' , which is the energy needed to excite 
respectively an electron into a higher-lying level, or a hole 
into a lower-lying level. For equidistant levels (S = S') 
these excitations are identical and we find n^a = np^. 

In Fig. [S] we show the same occupation numbers as in 
Fig. O but now calculated for different tunnel- couplings 
to the left lead. Again, the ground-state is partially 
depleted with increasing voltage, but due to the asym- 
metric couplings a clear asymmetry in bias is observed. 
For negative bias, only the p-states are being populated, 
due to a strong \s) — > \p) (cotunneling) amplitude pro- 
portional to tL2tR3 = 3. The \s) — > \h) amplitude is 
proportional to iLi^fi2 = 1/3 and therefore only hardly 
pumped at all. For positive bias, the situation is re- 
versed since (p|i?int|s) oc i^a^m = 1 is now much smaller 
than {h\Hint\s) oc tRitL2 — 3, whereby it is Uh which 
rises with voltage. Notice also that since {h\Hi^t\p) = 
there are no direct transitions between the two excited 
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states, whereby an incipient occupation of say, will 
be accompanied by a deletion of \s) and thereby further 
prohibit the pumping of \p). 

For the asymmetric tunnel-couplings chosen for Fig. [SI 
one can also observe a population inversion at positive 
bias. That is, for eV ^ S + g^sB, n^-^ increases rapidly 
and stays larger than the ground-state occupation ti^i 
for an extended bias-range. The positive bias demands 
electrons to predominantly jump from the left lead onto 
the dot and from the dot into the right lead. Therefore, 
^L2 ^ iii will make the \s) \h) much more likely than 
the \h) — > \s) transition at positive bias and the dot gets 
stuck in the excited state \h). As we shall argue later, 
this excited state does not sustain as high a current as the 
ground state and therefore this population inversion will 
in fact lead to regions of negative differential conductance 
(NDC). 

The 'direct' threshold for populating \h, f), say, is given 
by eV — 6. Nevertheless, from Fig. [5] it is seen that 
Uhi starts growing already at eV — S — g^sB. This 
slightly lower 'indirect' threshold is due to a cascade effect 
in which |s,J,) is populated for V > gfJ-sB and from 
there the threshold to t) is now only 5 — gnBB instead 
of 5. Further cascade effects are observed a.t eV = S 
(instead of the direct threshold S + g^isB) which is the 
threshold for the transition |s,i) \h,i). A similar 
effect for np^^ is observed at negative bias, with indirect 
thresholds already at eV — —{S'— gusB) and eV = 
—6'. Notice that these effects are just barely visible for 
symmetric couplings as in Fig. [51 As we shall see in the 
next subsection, these cascade effects show up as small 
steps in the nonlinear conductance. 



B. Cotunneling conductance 

In Fig. [7| we show the differential conductance for the 
same parameter set as in Fig. [5] As expected, we find 
a step in dl/dV whenever a new level enters the volt- 
age window. Furthermore, as pointed out in Refs. [l5ll20l . 
the voltage-dependence of the nonequilibrium occupation 
numbers promotes these conductance steps to pointed 
cusps. Since we are limiting our calculations to second 
order perturbation theory, these plots do not take into 
account the possible Kondo-enhancemcnt of these cusps. 
Nevertheless, the tendency is known from previous cal- 
culations: each cusp in the differential conductance to 
second order is logarithmically enhancecU^ and these log- 
arithmic divergences are contained roughly by the inverse 
life-time of the excited state involved in the relevant in- 
elastic cotunneling process2ii2^. 

1. Cascade induced side-peaks 

Even for the symmetric couplings in Fig. [71 we observe 
extra steps reflecting the aforementioned cascade effect. 
For example, the red (dashed) curve with gfisB = 0.25(5 
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FIG. 6: (color online) Occupation numbers of the spin, hole 
and particle states in the 3-leveI Anderson model for non- 
symmetric coupling to the leads, i.e. tL2/tRn = 3, tLs/tRn — 1 
and tLi/tiin = 1/3 where tun = G.l^U /vp. With the same 
finite magnetic field B/5 = 0.25, level splitting 5' /5 = 2 and 
temperature T/5 = 0.01 as in Fig. [5] we observe an occupa- 
tion inversion around eV ~ S + B since rihi; becomes higher 
occupied than the ground state Usi. This blocking of transi- 
tions and the cascade effect appearing at 5' ' — gfisB are dis- 
cussed in the main text. Gridlines refer to various thresholds: 
±gfiBB (blue), {5,5±gfiBB} (black), and {—S',—S'±gfiBB} 
(red). 



shows extra steps both at eV = S — g^sB (|s, |) — > t)) 
and at eV — 5' — g^sB (|s, J,) — > \p, t))- Notice that the 
cascade effect is only visible for B ^ 0, and disappears 
when 2gfiBB > 6(5') since n^j(npj) is then occupied 
before Ug^i. 

More investigations are necessary in clarifying to what 
extent these cascade-features in the conductance will be 
enhanced by the Kondo-effect, but even a small addi- 
tional step as seen in Fig. [3 can have important bear- 
ings for interpreting experiments in which the internal 
dot, or molecule states are not known in advance. Judg- 
ing from the magnetic field dependence, the three tran- 
sitions at eV = 6 — giiBB,6,6 + gusB seen in Fig. [71 
(imagine for a moment that S' ^ S) could in fact be mis- 
interpreted as direct transitions between a spin-doublet 
ground state and an excited S=3/2 state, as illustrated 
by the energy-diagrams in Fig. [51^^ This example shows 
that a spin-spectrum read off from B-dependent inelas- 
tic cotunneling-lines in a diamond-plot (dl/dV vs. Vg 
and Vsd) should be interpreted with some care, especially 
when enhancing the lines by plotting higher derivatives 
of the current. 



2. Negative differential conductance 

In Fig. [HI we plot the conductance for the parameters 
used in Fig. [SI and with varying values for the tunneling 
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FIG. 7: (color online) Dependence of the differential con- 
ductance dl/dV on the applied bias voltage in units of 
eV/5 for various strengths of the magnetic field gfisB/S = 
0, 0.25, 0.5, 0.75. The quantum dot is coupled symmetrically 
to the lead with the parameters given in Fig. [5] The kinks at 
±{5 — gfj,BB) and ±{S' ~ g^sB) originate from the occupation 
of excited states due to the cascade effect. 



amplitudes from the left lead to orbitals 1 and 2, respec- 
tively. As in Fig.ini we observe a pronounced asymmetry 
in bias-voltage, and with the largest difference in cou- 
plings to orbitals 1 and 2 (black/solid curve), we observe 
a sudden drop to negative differential conductance for 
eV > S. As mentioned before, the NDC is driven by the 
population inversion between \s) and \h). As argued in 
the previous section, the positive bias drives the dot into 
the excited state \h), and since (for these parameters) 
this state leads to a lower current than |s), the current 
decreases and we observe the NDC. The poor cotunneling 
across the dot in the state \h) is a simple consequence of 
the double occupancy of the well-coupled orbital 2, which 
demands that orbital 2 is emptied to the right lead be- 
fore it can be filled from the left, thus making use of the 
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FIG. 8: Comparison of direct and indirect transitions with 
a doublet groundstate and an excited doublet or quartet, re- 
spectively, at finite magnetic field. Cascade effect provides an 
indirect transition in the former, at energy S — gfisB, which 
matches a direct transition in the latter. 
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FIG. 9: (color online) Differential conductance dl jdV of the 
non-symmetric 3-level Anderson model with 3 electron for the 
parameters given in Fig. [B] The four different set of param- 
eters correspond to the ratio of tL2/tRn = {3, 2, 1.5, 1} with 
iLi Afln = l/(ti2 Afln) = {1/3, 1/2, 2/3, 1}, respectively. The 
parameter set with tL2/tRn ~ 3 and tLi/tHn = 1/3 corre- 
sponding to Fig. [6] shows NDC at the voltage values corre- 
sponding to an inversion of occupation probabilities. 

strong coupling = 3. In contrast, starting in \s) one 
can also fill in an electron in orbital 2 from the left lead 
first and then tunnel out into the right lead from either 
orbital 2 or from orbital 1, thus allowing more current- 
carrying tunneling processes involving 1^2 = 3. 

Thus a sufficiently strong difference in coupling to two 
of the orbitals will cause a population inversion with a 
concomitant decrease in the current as bias is increased. 
For intermediate coupling asymmetry, we see from the 
other curves in Fig. [5] that the NDC disappears while an 
asymmetry in bias remains. 



IV. DISCUSSION OF EXPERIMENTS 

A. Experiments on Carbon Nanotube 

The data shown in Fig. [3] were recorded on the same 
sample as was discussed in Ref. M. In fact, this is the 
neighboring charge-state to the even number occupied 
state studied in that work. In Ref. iS, some of the present 
authors obtained a very gratifying fit to the inelastic co- 
tunneling line reflecting a strong transition from singlet 
ground- to an excited triplet state. In particular, it was 
argued that the very sharp cotunneling line could only 
arise from the joint effect of nonequilibrium pumping 
(finite-bias occupations) of the triplet state and substan- 
tial logarithmic enhancements from the nonequilibrium 
Kondo-effect. 

As already stated, we shall not embark on higher 
order perturbation theory in this paper, which means 
that we should not expect to obtain any quantita- 
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FIG. 10: Modeling ol the CNT QD analyzed in Ref. i 
as measured in the neighboring diamond, i.e. Fig. O 
where orbital scattering is observed at 5 = 1.5 meV. 
Further parameters of the plot are {t^i, tL2, iiii, tfl2} ~ 
{0.032, 0.028, 0.108, 0.063} ^17/ ;/f, B = 0, 0.14, ... 1.4 meV, 
and T — 81mK ~ 0.008meV. The conductance is given in 
arbitrary units and the curves are offset for clarity. 



FIG. 11: (a) and (b) Plot of dl /dV vs. gate, and bias-voltage, 
for the InAs-wire quantum dot with _B = T and 0.25 T, 
respectively. Strong tunnel-coupling perturb the Coulomb- 
blockade diamond and particularly strong sequential tunnel- 
ing lines are clearly present (thick dark lines) in the region 
with a mixture of 2 and 3 electrons on the dot. c) Calcula- 
tion of the sequential tunneling lines from solving semiclas- 
sical rate-equations (conductance in arbitrary units) bring- 
ing out some of the main lines and indicating the states 
which are involved in the dominant transitions and which 
are involved in the cotunneling for A*' = 3. Parameters are 
5 = 1.4 meV, 5' = 1.1 meV, T = 300mK, U = Z meV, 
thi = tm = tL3 = tR3 = 1, tL2 = 0.4, and tR2 = 2.5. 



tive agreement with the data in Fig. [3] Nevertheless, 
since basically all parameters were fixed by the fit in 
Ref. H, we shall assume that these remain largely un- 
changed when passing on to the neighboring charge-state 
and use them as input for a numerical evaluation of 
the nonlinear conductance there. The result is shown 
in Fig. \W\ and the gross features such as bias-voltage 
asymmetry and B-field splittings match those of Fig. [3] 
quite well, albeit with a complete lack of sharp finite- 
bias peaks, as expected. The parameters from Ref. i8| 
were tunnel couplings of the two lowest lying levels of 
{tLi,tL2,tm,tR2} = {0.032, 0.028, 0.108, 0.063} ^tZ/i^F 
(see Fig. 4 in Ref. 8), an orbital splitting oi S « 1.5 
meV, a charging energy of U « 3.0 meV, and a higher 
lying level at(5' = A — (5~3.1 meV, A being the single- 
particle level spacing in the nanotube. Since S' > U, 
charge-excitations set in before orbital 3 is reached and 
this latter orbital is therefore not included in the calcu- 
lation. 



B. Experiments on InAs-nanowires 

The data shown in Fig. [5] constitute the first observa- 
tion of Kondo-effect in InAs-nanowire quantum dots and 
were discussed by some of the present authors in Ref. [^. 
Already in that paper the 3-orbital model was invoked to 
explain the finite-bias peaks flanking the zero-bias Kondo 
peak and, as mentioned in the introduction, this exper- 
iment was the main motivation for this more detailed 
investigation of that model. 

Apart from the zero-bias Kondo peak, it is the sharp 
finite-bias peaks which dominate these data: A very 



strong peak close to eV = —1.4 meV together with 
a somewhat weaker peak close to eV = 1.1 meV 
(cf. Fig. [2]). This asymmetry in bias- voltage was not 
addressed in Ref. d, but as we shall argue it can read- 
ily be understood in terms of an asymmetry in tunneling 
amplitudes of orbital 2 to source and drain. The excited 
state giving the strong sequential tunneling line in the 
upper left corner of Fig. [TT] (a) is ascribed to the excited 
3-particle state with one electron in orbital 3 instead of 
orbital 2. This line connects to the inelastic cotunneling 
line at positive bias inside the N=3 diamond (|s) — > \p) 
transition). From this we infer that S' w 1.1 meV. No- 
tice that a sequential tunneling line corresponding to a 
transition from the N=2 ground-state to the N=:3 \h)- 
state is only possible to higher order in the tunneling 
amplitudes and is therefore strongly suppressed in the 
experimental data. The inelastic cotunneling line at neg- 
ative bias can now be ascribed to an (N=3) \s) — * \h) 
transition at energy 5 ~ 1.4 meV. This is again consis- 
tent with the lower left corner of Fig. [TT] (a) showing that 
this line connects to a very strong sequential tunneling 
line which must correspond to a transition from the N=4 
ground-state to the N=3 |/i)-state. At this point it is 
the (N=3) |p)-state which is not directly coupled to the 
N=4 ground-state and hence suppressed in the data. In 
Fig.[TT](c), we show the result of solving a set of semiclas- 
sical rate-equations. The dark lines correspond to high 
conductance (arbitrary units) due to sequential tunnel- 
ing only. We do not attempt a 6 parameter fit of Fig. [TT] 
(a), and tunneling amplitudes are therefore chosen with 
a very simple asymmetry, assuming that only orbital 2 
is coupled in a special way. This choice of parameters 
captures some, but not all, of the gross features of the 
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FIG. 12: (color online) Plot of dl/dV vs. bias voltage, 
V , calculated within second order perturbation theory. Pa- 
rameters are chosen so as to approximately match the volt- 
age asymmetry brought out by the experiment (compare 
Fig. (TT)). Parameters are B = 0, 5 = 1.4 meV, 5' = 1.1 
meV and T = 0.001 meV (blue/solid curve). The black 
(dashed-dotted) curve corresponds to all equal tunneling am- 
plitudes and the red (dashed) curve is the same as the blue 
(solid) but with temperature corresponding to the experiment 
(Toxp = 300mK~ 0.025meV). 



FIG. 13: (color online) Plot of dl/dV vs. bias voltage, V, 
calculated within second order perturbation theory for sym- 
metric coupling to the leads and neglecting the particle state 
{5' > Ec). The g-factors are different for the two involved 
levels: g2 = 2gi. The orbital peak moves towards higher 
voltages such as ~ (1 — (71/(72 )g2/^sS/2 with increasing mag- 
netic field g2tJ^BB/5 = {0,0.25,0.75}. The red (dashed) 
curve at g2fiBB/5 = 0.75 illustrates that the shift is also 
observable at a temperature corresponding to the experiment 
(Texp = 300mK~ 0.025 meV). 



data in Fig. [TT] (a). 

The question remains why there is only one and not 
two inelastic cotunneling peaks at both positive and neg- 
ative bias. As demonstrated by the calculation shown in 
Fig. 1121 however, the observed asymmetry in bias- voltage 
can be understood quite simply as orbital 2 being asym- 
metrically coupled to the source and drain electrodes. 
Notice that it is the nearly equidistant levels in the InAs- 
wire quantum dot which prompts us to incorporate the 
effects of both h and p excited states. For asymmetric 
couplings the 6' — l.lmeV and S — 1.4meV peaks in 
Fig. [12] dominate for positive and negative voltages, re- 
spectively, in contrast to symmetric coupling where both 
can be observed. The small steps observed for zero tem- 
peratue (blue/solid curve) are washed out for large tem- 
perature smearing (red/dashed curve) and the nonlinear 
conductance is now asymmetric with respect to voltage. 
For a larger difference in the effective 6 and 6' , as is typi- 
cal for carbon nanotube quantum dots where d would be 
the subband-splitting and 6' would be of the order of the 
single-particle level-spacing (A = 6' + 6), the \p) state 
would most often be masked by charge-fluctuations, i.e. 
S' will often be comparable to the charging-energy, Ec- 

As demonstrated recently in the experiment by Csonka 
et alr^, the strong spin-orbit coupling in InAs can give 
rise to very different g-factors for neighboring quantum 
dot orbitals or energy-levels. In that experiment this 
was brought out particularly clear by a modulation of 
the orbitals by an additional top-gate. Fixing the top- 
gate voltage and adjusting the back-gate, it was shown 



that two neighboring levels could have g-factors of 1.9 
and 10, respectively. It is interesting to note that such a 
difference in g-factors would give rise to an apparent in- 
dependence of the level-spacings. In our model, we would 
expect the following _B-dependent inelastic cotunneling 
thresholds: 

6{B) = Eh^ -Es^=S+ (.ga - gi)fiBB/2, 
5'{B) = Ep^ - Es^- =6' + {g2 - 53)/^sB/2. 

In Ref. ,5i, a detailed investigation of the splitting of both 
the zero-bias Kondo peak and the finite(positive)-bias co- 
tunneling peak with magnetic field indeed revealed two 
different g-factors. Interpreting the finite-bias cotunnel- 
ing as above, we can thus infer from Ref. 5 that orbital 2 
has (72 = 7.7 and orbital 3 has (73 = 8.5. This difference 
is rather small and would cause a largely negligible shift 
of roughly 0.02 meV of 5' at largest applied fields (0.9 
Tesla). In Fig. [l3] we show the result of a calculation 
with slightly different g-factors by a factor of 2. This is 
the cause of the slight shift of the main kink in the blue 
(solid) curve to a value above eV = 5. The observation of 
such i?-dependent cotunneling thresholds thus provides 
an interesting consistency check on the difference in g- 
factors of two neighboring orbitals. 

In closing, we note that the diamond-plot in the upper 
left panel of Fig. [TT] has a few strong-coupling irregu- 
larities which we do not address in this paper. First of 
all, the inelastic cotunneling lines display a slight depen- 
dence on gate-voltage. This is a feature which was in- 
vestigated in detail for a carbon nanotube quantum dot 
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in Ref. '9', where it was explained in terms of second or- 
der tunneling-renormalization with energy denominators 
depending on gate-voltage. Second, the inelastic cotun- 
neling line at positive bias, which we assigned to the ex- 
cited p-state, appears to continue straight through the 
sequential tunneling region mixing N=3 and N=4 states. 
At this point we merely speculate that this is related to 
a very strongly coupled orbital 3, which might allow the 
partially populated N=4 states to conduct via virtual 
transitions to the N=5 state. How, and if, this would 
work out in detail will be assessed in future work. 



V. CONCLUSIONS 

In this paper we have analyzed the inelastic cotunnel- 
ing spectroscopy of a system with three electrons dis- 
tributed among three orbitals. Starting from the 3- 
orbital Anderson model in the Coulomb blockade regime, 
we have projected out charge-fluctuations and derived an 
effective cotunneling (Kondo) model describing elastic as 
well as inelastic tunneling induced transitions between 
different orbital and spin states of the impurity or quan- 
tum dot in the N=3 charge-state. By treating all tunnel 
couplings independently we can study asymmetric behav- 
ior inherent to most of the experimental samples. 

All calculations were performed within second order 
nonequilibrium perturbation theory, and the bias depen- 
dent nonequilibrium occupations of the impurity-states 
were shown to give rise to marked cusps at voltages 
matching the different transitions. We demonstrate how 
a bias induced cascade in the occupations can lead to 
indirect transitions between two excited states. Further- 
more, we find that certain asymmetries in the tunneling 
amplitudes of the three orbitals can lead to a population 
inversion, which in turn can be accompanied by lines of 
negative differential (cotunneling) conductance (NDC). 

We have revisited the carbon nanotube data from 
Ref. m with focus on an odd occupied charge-state. Im- 
porting the tunnel-couplings inferred from Ref. we 
confirm the asymmetry in the conductance strength at 
positive and negative bias. These parameters lead to 
steps at finite-bias at the orbital transitions, and we be- 
lieve that the experimentally observed peaks are mainly 
due to logarithmic enhancements from a nonequilibrium 
Kondo-effect. 

We have also analyzed the InAs-wire data from Ref. 
in greater detail and provided an explanation of the bias- 
asymmetry in peak-positions. This was understood as a 
slight difference in the two level- spacings together with an 
asymmetry in the tunnel-coupling of orbital 2 to source 
and drain electrodes. This asymmetry in couplings was 
roughly consistent with the overall stability diagram for 
the N=3 diamond, as confirmed by a semiclassical rate- 
equation calculation. Finally, we have pointed to an ad- 
ditional consequence of strong spin-orbit coupling in the 
InAs-wires, namely the apparent -B-dependence of the 
level-spacing caused by different g-factors for different 



orbitals. In the present experiment this difference was 
too small to give a noticeable effect, but with differences 
as quoted in Ref. [l^ it should be a pronounced effect^i. 
Possible further influence of strong spin-orbit coupling on 
the cotunneling spectroscopy of inelastic spin and orbital 
transitions constitutes an interesting problem on its own 
right and will be addressed in a separate publication. 
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APPENDIX A: SCHRIEFFER- WOLFF 
TRANSFORMATION AND COUPLING 
CONSTANTS 

In this work we concentrate on the physics inside a 
Coulomb diamond where the 3-level quantum dot is oc- 
cupied by exactly three electrons. All virtual processes 
to two or four electron states can be treated perturba- 
tively and lead to a Kondo-like interaction and potential 
scattering terms. 

Using a Schrieffer- Wolff transformation^'^ we find to 
second order in the tunneling Hamiltonian Htun the in- 
teraction Hamiltonian 

^int = ^^amtpn'^aafmal^ _ ^ fna'^^a' 



cr,cr 



Eaa — -E'4e- 



which can be rewritten as 



aml'0n ''fSn'' 



+ 3U em + 2U 



a. a' 



e„ + 3C/''""'''5'^' 



(Al) 



The Kondo exchange interaction for the levels n,m = 
{1, 2, 3} and thus for the orbital states a,b = {s, h,p} can 
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straightforwardly be read off this expression 

Til _ jhh _ 9/ 

<: 



tion Hamiltonian 



t22 _ TSS _ 91 



t33 _ tPP _ 01 



t12 _ TS/l 



jZL _ jhs 
— •-'aP 



t32 _ TP'' — of 



t23 _ Tsp _ 
— -Jap — 





i pi^al 


El + 3U 


Si + 2U 


ta2tp2 




62 + 3U 


S2 + 2U 


ta3t*p3 


*/33*"3 


£3 + 3t/ 


£3 + 2C/ 




tpitci2 


El + 3U 


62 + 2U 


tal't*p2 


t*p2^al 


e2 + iU 


ei+2U. 


ta2t*pz 


t*p3^a2 


£3 + 3C/ 


e2 + 2U 


ta3t*p2 




e2 + 3U 


e3 + 2U 



which leads to the interaction Hamiltonian as defined in 
the main text. Note that the potential scattering which 
changes the orbital index is of the same strength as the 

spin Kondo scattering with an orbital change, j'^^p''^- In 
lowest order, i.e. involving two hopping processes, there 
is no transition between the hole and the particle state 
and we assume for the following 



jhp _ jph _ 
•^aP — '^ap — 



Furthermore we disregard a constant shift due to the 
following potential scattering contribution 



pot. scat 



yCpc 



ap 



since it is negligible small at the particle-hole symmetric 
point £2 = —5/2 U, 



Gap — ( — 

V£i 



f3C/ £i + 2C/ 

ta2tfj2 tp2't 



2UJ V 



+ 



.£2 + 3C/ e2 + 2UJ \e3 + 3U £3 + 2C/ 

-^allpi _ J2 'a3l/33 



■(f//2)2 - {b'f 



This constant offset in the current of order J/f/^ is ne- 
glected in our calculation. 



APPENDIX B: CALCULATION OF THE 
SELF ENERGIES 



The Green's functions for the orbital states a,b = 
{s,h,p} are calculated by an expansion in the interac- 



9ab,aa'{r,r') = -i (Tc,, d,. (r)^, (t' ) )^ 

Jck ^ ' 

^"'^ I' dTidT2(Tc,.d,,(r)i?i„t(Ti)i/i„t(r2)dL,(r') 



Ck 



(Bl) 



The linear order does not contribute since it contains the 
expectation value of (sap) which is zero if the spin in 
the leads are not polarized. In the case of ferromagnetic 
leads this order has to be taken into account, but in our 
setup the leading order is second order in the coupling. 

We get contributions from the spin part of the inter- 
action Hamiltonian and of the potential scattering part. 
A mixing between both parts does not appear in second 
order due to the convolution of the leads contribution. 

For example the spin part is given by (Einstein's sum 
rule) 

g(2),Spin. M f I ^, 

X (^CkcLj i^l)Cpa^ (n)cl,,, {^2)0/^,^^ {T2)) ^ 

x(Tc,d,,(r)J-''^Sl;,^(n)J:?|?S-,Jr2)dL,(r'; 
The conduction electron spins contract to 

^(0 
^Pc 



^a'aTaa'6a'pSaP'G''^l,{Ti,T2)G'-°}{T2,Ti). 



Since the leads are not magnetic, the conduction electron 
Green's functions does not depend on spin and the sum 
over the spins acts only on the r-matrices and we can 
thus introduce the conduction electron spin susceptibility 

XapiT2,Tl) = Ga^(r2,Ti)G^°^(Ti,r2). 

As discussed before there can in general be solutions 
with a ^ b which we neglect in this setup>i^. The re- 
sult for the diagonal Greens 's function yields finally the 
second order self energy by comparison with the Dyson 
series: 

^aa{Tl,T2) = jQ^'^'^pJpa + -^q^J-^Iq) Xa/3 (t2 , n ) 
[2Xcrcr'Tcr'CT + ^cr,a''?a.b]Gba' {ti , T2) 

+ ^t{ri,T2) 

where the inelastic and elastic potential scattering con- 
tributions are defined as 

^ab = {Sa,sSb,p + Sb,sSa,p + Sa,sSb,h + Sb,sSa,h) 
1 



Kairi,r2) 



8^ 



jaa jaa^, ^ \ 

JapJpaXaP[T2,Ti) 



Note that in the rate equation E^'^ does not contribute 
since it does not contain transitions between states. 
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The self-energy components needed in Eq. (|19p are now Langreth rules, 
readily obtained from analytical continuation using the 
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